DtO  251975 


AEDC-TR-75-118 


.2 


CALCULATION  OF  THE  BOUNOARY-LAYER 
GROWTH  IN  A  LUDWIEG  TUBE 


VON  KARMAN  GAS  DYNAMICS  FACILITY 
ARNOLD  ENGINEERING  DEVELOPMENT  CENTER 
AIR  FORCE  SYSTEMS  COMMAND 
ARNOLD  AIR  FORCE  STATION,  TENNESSEE  37389 


December  1975 


Final  Report  for  Period  July  1973  —  December  1974 


Approved  for  public  release;  distribution  unlimited. 


i-w-  iijf  o£  'J-  ^ 

F4CG00-7D -0-0001 


Prepared  for 

DIRECTORATE  OF  TECHNOLOGY  (DY) 

ARNOLD  ENGINEERING  DEVELOPMENT  CENTER 
ARNOLD  AIR  FORCE  STATION,  TENNESSEE  37389 


4. 


NOTICES 


When  U.  S.  Government  drawings  specifications,  or  other  data  are  used  for 
any  purpose  other  than  a  definitely  related  Government  procurement 
operation,  the  Government  thereby  incurs  no  responsibility  nor  any 
obligation  whatsoever,  and  the  fact  that  the  Government  may  have 
formulated,  furnished,  or  in  any  way  supplied  the  said  drawings, 
specifications,  or  other  data,  is  not  to  be  regarded  by  implication  or 
otherwise,  or  in  any  manner  licensing  the  holder  or  any  other  person  or 
corporation,  or  conveying  any  rights  or  permission  to  manufacture,  use,  or 
sell  any  patented  invention  that  may  in  any  way  be  related  thereto. 


Qualifled  users  may  obtain  copies  of  this  report  from  the  Defense 
Documentation  Center. 


References  to  named  commercial  products  in  this  report  are  not  to  be 
considered  in  any  sense  as  an  endorsement  of  the  product  by  the  United 
States  Air  Force  or  the  Government. 


This  report  has  been  reviewed  by  the  Information  Office  (01)  and  is  releasable 
to  the  National  Technical  Information  Service  (NTIS).  At  NTIS,  it  will  be 
available  to  the  general  public,  including  foreign  nations. 


APPROVAL  STATEMENT 


This  technical  report  has  been  reviewed  and  is  approved  for  publication. 


FOR  THE  COMMANDER 


CARLOS  TIRRES 
Captain,  USAF 
Research  &  Development 
Division 

Directorate  of  Technology 


c? 

ROBERT  O.  DIETZ 
Director  of  Technology 


UNCLASSIFIED 


REPORT  DOCUMENTATION  PAGE 

READ  INSTRUCTIONS 

BEFORE  COMPLETING  FORM 

1  report  number  2.  GOVT  ACCESSION  NO. 

AEDC-TR-75-118 

3.  RECIPIENT’S  CATALOG  HUMBER 

4  title  fand  Subr/tfe^ 

CALCULATION  OF  THE  BOUNDARY-LAYER  GROWTH 

IN  A  LUDWIEG  TUBE 

S  TYPE  OF  REPORT  6  PERIOD  COVERED 

Final  Report-July  1973  - 
December  1974 

6.  PERFORMING  ORG.  REPORT  NUMBER 

^  AuTHORfsJ 

James  C.  Sivells,  ARO,  Inc. 

8.  CONTRACT  OR  GRANT  NUMBERfa; 

1 

9  PERFORMING  organization  NAME  AND  ADDRESS 

Arnold  Engineering  Development  Center  (DY) 
Arnold  Air  Force  Station,  Tennessee  37389 

10  PROGRAM  element.  PROJECT.  TASK 
AREA  A  WORK  UNIT  NUMBERS 

Program  Element  65807F 

n.  controlling  office  name  and  address 

Arnold  Engineering  Development 

Center  (DYFS) 

Arnold  Air  Force  Station.  Tennessee  37389 

1Z.  REPORT  DATE 

December  1975 

•  3.  NUMBER  OF  paces 

36 

14  monitoring  agency  NAME  ft  ADDRESSf/l  di/f0rBnl  from  Contfaltlng  OHicb) 

IS  security  class,  (el  ihle  report) 

UNCLASSIFIED 

1  1  |i|  1  H  II  1 

16  distribution  STATEMENT  ref 

Approved  for  public  release;  distribution  unlimited. 

1 

17  DISTRIBUTION  STATEMENT  (ol  Me  mbttrmcl  mnfnd  In  Block  20,  11  dlllortnl  Irom  Ropart) 

16  SUPPLEMENTARY  NOTES 

Available  in  DDC 

19  KEY  WORDS  (Contlnum  on  raveraa  alda  It  naceaaaty  and  identity  by  block  nutnber) 

computations  Ludwieg  tube  expansion 

boundary  layer  transonic  flow  wave 

growth  wind  tunnels  (pilot) 

measurement  Reynolds  number 

20  A  BSTRACT  (Continue  on  revmrae  eide  It  neceoeery  and  Identity  by  block  number) 

Experimental  boundary-layer  measurements  obtained  in  a  Ludwieg 
tube  used  to  drive  a  pilot  transonic  tunnel  were  compared  with  . 
values  calculated  by  a  procedure  developed  by  E.  Becker  and  with 
those  calculated  by  a  method  containing  several  modifications  to 
Becker's  method.  The  modifications  fall  into  three  general 
categories:  the  use  of  a  skin-friction  law  and  velocity  profile 
exponent  which  are  more  accurate  at  high  Reynolds  numbers; 

DD  I  JAN *73  1473  edition  OF  I  NOV  65  IS  OBSOLETE 


UNCLASSIFIED 


UNCLASSIFIED 


20.  ABSTRACT  (Continued) 


treatment  of  the  momentum  equation  as  axisymmetrlc  instead  of 
two-dimensional;  and  calculation  at  a  specified  location  other 
than  at  the  origin  of  the  centered  expansion  wave.  Inasmuch  as 
these  modifications  greatly  Improved  the  agreement  with 
experimental  values,  they  are  presented  herein. 


UNCLASSIFIED 


AEDC-TR-7S-118 


PREFACE 

The  work  reported  herein  was  conducted  by  the  Arnold  Engineering  Development 
Center  (AEDC),  Air  Force  Systems  Command  (AFSC),  under  Program  Element  65807F. 
This  work  was  done  by  ARO,  Inc.  (a  subsidiary  of  Sverdrup  &  Parcel  and  Associates, 
Inc.),  contract  operator  of  AEDC,  AFSC,  Arnold  Air  Force  Station,  Tennessee.  The 
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1.0  INTRODUCTION 

In  a  wind  tunnel  driven  by  a  Ludwieg  tube  (Ref.  I),  the  air  in  the  wind  tunnel 
and  tube  is  initially  compressed  to  a  desired  charge  pressure.  Flow  is  initiated  by  breaking 
a  diaphragm,  or  quickly  opening  a  valve,  located  downstream  of  the  test  section.  As  the 
air  is  released,  an  expansion  wave  is  created  and  travels  upstream  to  the  closed  end  of 
the  tube  where  it  is  reflected  and  returned  to  the  contraction  section  at  the  downstream 
end.  Due  to  viscous  effects  in  the  airflow  generated  by  the  expansion  wave,  a  boundary 
layer  is  formed  whose  thickness  increases  with  time.  During  this  excursion  of  the  expansion 
wave,  the  stagnation  pressure  of  the  central  core  of  flow  through  the  wind  tunnel  is 
essentially  constant  until  the  thickness  of  the  boundary  layer  at  the  downstream  end  of 
the  tube  approaches  the  radius  of  the  tube.  Thus,  the  useful  run  time  for  the  Ludwieg 
tube  ^^d  tunnel  depends  upon  the  initial  air  temperature  which  determines  the  velocity 
of  the  head  of  the  expansion  wave,  the  length  of  the  tube  which  is  traversed  by  the 
expansion  wave,  and  the  diameter  of  the  tube  which  determines  the  velocity  of  the  air¬ 
flow  and  relative  to  which  the  boundary-layer  thickness  eventually  becomes  critical. 

Shortly  after  the  conception  of  the  tube  wind  tunnel,  calculations  of  the  growth 
of  the  boundary  layer  were  made  by  E.  Becker,  (Refs.  2  and  3).  Values  calculated  by 
Becker's  method,  however,  considerably  underestimated  the  boundary-layer  thicknesses 
obtained  experimentally  in  a  pilot  tunnel  as  shown  in  Fig.  1,  which  was  presented  in 


AEDC-VKF  I  mproved  Theory 
Becker's  Theory 


Figure  1.  Charge  tube  boundary-layer  thickness. 
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Ref.  4.  In  this  figure,  the  Reynolds  number  is  based  on  the  distance  Ax  between  the 
head  of  the  expansion  wave  and  the  measuring  station.  Data  are  presented  for  two  times, 
when  the  length  of  the  expanaon  wave  AW  (head  to  tail)  was  about  SS  and  also  about 
75  percent  of  Ax.  Three  values  of  charge  pressure  P4  are  indicated. 

Scrutiny  of  Becker's  procedure  indicated  three  categories  in  which  modifications  could 
be  made  to  improve  its  correlation  with  the  experimental  data;  the  use  of  a  skin-friction 
law  and  boundary-layer  velocity-distribution  which  are  more  accurate  at  high  Reynolds 
number;  treatment  of  the  momentum  equation  as  axisymmetric  instead  of  two-dimensional; 
and  calculation  at  a  specified  location  other  than  at  the  origin  of  the  centered  expansion 
wave.  The  latter  is  particularly  important  for  a  transonic  wind  tunnel  utilizing  a  plenum 
chamber  surrounding  a  porous-wall  test  section  because  the  establishment  of  steady-state 
conditions  in  the  test  section  determines  when  the  tail  of  the  expansion  wave  starts 
upstream  through  the  Ludwieg  tube. 

When  Becker's  procedure  was  so  modified,  the  results  indicated  as  the  AEDC-VKF 
improved  method  in  Fig.  1  were  obtained.  The  improvement  in  correlation  with  data  from 
the  pilot  tunnel  was  considered  to  be  sufficient  to  warrant  the  use  of  the  modified  method 
in  any  future  applications. 


2.0  THEORY  OF  OPERATION 


The  principle  of  operation  can  be  described  with  the  help  of  the  wave  diagram  of 
Fig.  2.  This  diagram  differs  from  the  usual  diagram  which  considers  the  diaphragm  to 
be  located  at  the  downstream  end  of  the  tube.  In  the  practical  case  of  a  transonic  wind 
tunnel,  the  diaphragm  or  start  valve  is  located  downstream  of  the  test  section.  When  flow 
is  initiated,  the  expansion  wave  travels  upstream.  The  velocity  of  the  head  of  the  wave 
is  the  speed  of  sound  ao  at  the  temperature  of  the  charge  air.  The  velocity  of  the  air 
leaving  the  tube  uj  is  determined  by  the  ratio  of  the  tube  area  to  the  area  where  the 
flow  is  sonic. 


Atube 

A* 


j_  /_!_  Ill 

\7+  1  ^  7+1 


T-t-l 

2(7-1) 


(1) 


and 


Ml  =  ui/ai  (2) 

The  velocity  of  sound,  ai ,  in  the  air,  after  the  flow  is  established,  is  related  to  ao. 


(3) 


6 


AEDC-TR-75-118 


The  velocity  of  the  tail  of  the  wave  is  ai  -  ui ;  thus,  the  length  of  the  expansion 
wave  increases  with  time.  The  wave  length  is  approximately  linear  with  time  after  the 
Avave  is  completely  in  the  charge  tube.  For  simplicity,  it  is  also  shown  in  Fig.  2  as  linear 
downstream  of  the  charge  tube  by  extending  the  lines  denoting  the  head  and  the  tail 
of  the  wave  until  they  meet  at  an  effective  location  of  the  origin  of  the  expansion  wave. 
In  the  actual  case  of  a  porous- wall  transonic  test  section,  some  time  is  spent  in  establishing 
steady-state  conditions  in  the  test  section  and  plenum  chamber.  This  in  turn  produces 
a  noncentered  wave  process  in  the  charge  tube.  The  start  time  indicated  in  Fig.  2  is  defined 
as  the  period  of  time  during  which  the  pressure  at  the  end  of  the  charge  tube  decreases 
from  its  charge  value  P4  to  the  value  of  Pi  following  the  passage  of  the  tail  of  wave, 
where 


27 

P1/P4  =  (ai/ao):^TT"  (4) 

From  a  knowledge  of  the  start  time,  either  from  measurement  or  estimate,  the  effective 
time  and  position  at  which  the  wave  length  is  zero  can  be  determined. 

In  the  development  of  the  boundary-layer  calculation,  Becker  solved  the  momentum 
equation  twice:  once  for  the  growth  of  the  boundary  layer  within  the  expansion  wave 
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and  again  for  the  growth  of  the  boundary  layer  behind  a  fictitious  concentrated  wave 
of  zero  length  moving  at  a  velocity  V  intermediate  between  the  velocities  at  the  head 
and  tail  of  the  actual  expansion  wave.  The  value  of  this  velocity  was  found  by  finding 
the  distance  behind  the  fictitious  wave  at  which  the  boundary-layer  thickness  was  the 
same  as  that  at  the  tail  of  the  expansion  wave  from  the  first  solution.  Becker  introduced 
the  variable  77,  where 


T?  = 


(5) 


and  17  varies  from  0  at  the  head  of  the  .expansion  wave  to  771  at  the  tail  and 

7+1  //  7-1  \ 

171  =  —  Mi/(l  +  —  Mij  (6) 

When  Becker  matched  his  boundary-layer  solutions,  he  obtained  a  relation  between  V 
and  ui  (or  ao).  When  his  relationship  is  expressed  in  a  series  on  powers  of  771,  it  can 
be  shown  that 


V/ao  =  1  -  2771 /3  (7) 

is  sufficiently  accurate  for  practical  values  of  771  up  to  about  0.45  and  is  independent 
of  the  choice  of  velocity  profile  exponent  and  friction  coefficient  law  since  these  are 
contained  only  in  the  negligibly  small  coefficients  of  771 2  and  higher  powers  of  771 .  In 
the  subsequent  development  of  the  modifications  to  improve  the  correlation  of  theory 
with  experiment,  it  is  assumed  that  this  relationship,  Eq.  (7),  is  accurate  inasmuch  as 
only  the  flow  behind  the  fictitious  concentrated  wave  is  considered. 


3.0  MOMENTUM  EQUATION 


The  time-dependent  momentum  equations  for  internal  tube  flow  is  adapted  from 
Ref.  5  as 


/  6  -  f")  (P  -  Pi)  dy  +  -^  (P1U161)  + 

o'  ' 


3ui  9 

P1U161  +  (Piui20i)  =  Tw  (8) 


•■■(•(.■■fX'-TSr)* 


(9) 

(10) 
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6i  =  6*  -6*^/2r  (11) 

di  =  d  -  9^I2t  (12) 

The  quantities,  5i  and  di,  may  be  considered  to  be  the  displacement  and  momentum 
thicknesses  when  the  boundary-layer  thickness  is  small  with  respect  to  the  radius  of  the 
tube.  When  the  boundary-layer  thickness  is  not  small  relative  to  the  radius,  the  true  values 
of  5*  and  d,  obtained  from  mass-defect  and  momentum-defect  considerations  must  be 
used  to  assess  their  effects  on  the  flow  but  5i  and  di  are  still  used  in  Eq.  (8).  Becker 
did  not  use  the  (I  -  y/r)  term  in  Eq.  (8),  (9),  or  (10),  in  which  case  Eq.  (8)  reduces 
to  the  two-dimensional  form.  In  Ref.  1,  Becker  omitted  the  first  term  of  Eq.  (8).  In 
Ref.  2,  he  included  this  term  and  also  included  the  effects  of  compressibility  and  heat 
transfer  on  the  ratio  of  p/pi  within  the  boundary  layer  and  upon  the  friction  coefficient, 
and  obtained  essentially  the  same  result  (within  about  one  percent)  as  for  the 
incompressible  case,  at  least  for  the  usual  conditions  of  operation  of  a  Ludwieg  tube. 

The  momentum  equation  can  be  simplified  if  it  is  assumed  that  the  free-stream  values 
of  ui  and  pi  are  independent  of  time  and  a  definition  is  added 


0 

Also,  in  Eq.  (8)  the  direction  of  x  is  positive  in  the  direction  uj ,  but  for  the  present 
purpose  X  is  positive  in  the  direction  of  wave  propagation.  Then 


1  a  bdi 

- (p*  +  6i)  - 

ui  at  ^ 


Tw 


Cf 

2 


ax  Piui^ 

It  is  assumed  that  the  ratio  (p*  +  6i)/0i  is  relatively  independent  of  time.  Then 


(14) 


p*  +  6i  a(9i  a^i  Cf 

^  —  (15) 

01  Uf  at  ax  2 

This  equation  can  be  integrated  through  the  introduction  of  the  distance  X  defined  as 

E  =  Vt  -  X  (16) 

which  is  the  distance  from  the  concentrated  wave  to  the  point  x,  and  x  and  t  are  zero 
at  the  effective  origin  of  the  expansion  wave.  It  is  further  assumed  that  the  friction 
coefficient  is  a  function  of  equivalent  flat-plate  momentum  thickness,  0c,  such  that 

Cf/2  =  d0c/dX  (17) 
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In  the  absence  of  a  longitudinal  pressure  gradient,  Eq.  (17)  can  be  integrated  to  give 

Cf/2  =  0c/X  .  (18) 

where  the  constant  of  integration  is  neglected.  Eq.  (15)  can  thus  be  integrated  to  give 


/  V  p*+5i  \ 

— "  V 


01  =  0c  =  X  Cf/2  =  (Vt-X)  Cf/2 


(19) 


After  rearranging 


(Vt.x)CF/2 

l+JL  elm-  (20) 

“1  fll' 


Finally, 


-  JL  (Vt  -  X) 

Bi  ^  ,  V  (21) 

“1  »! 

The  ratios,  6/0i ,  p*/0i,  6i/0i,  are  obtained  from  Eqs.  (9),  (10),  and  (13)  after  the  velocity 
and  density  distributions  are  assumed  and 

V/ui  =  (1/Mi)  -  (5  -  y)l6  (22) 

from  Eqs.  (2),  (3),  (6)  and  (7).  The  friction  coefficient  must  be  in  a  form  which  can 
be  integrated  with  respect -to  the  momentum  thickness  and  must  be  corrected  for  at  least 
first-order  effects  of  compressibility  and  heat  transfer.  The  heat  transfer  is  from  the  wall 
to  the  air  inasmuch  as  the  mass  of  the  wall  is  considered  to  be  sufficient  that  its  temperature 
is  essentially  constant  during  the  short  run  time. 

Except  for  the  definition  of  the  various  parameters,  the  above  derivation  follows 
that  of  Becker.  When  x  =  Vt,  i.e.  at  the  location  of  the  concentrated  wave,  the 
boundary-layer  thickness  is  zero  and  increases,  for  a  given  value  of  t,  as  x  decreases  to 
zero.  In  the  present  treatment,  the  lowest  value  of  t  to  be  considered  is  the  time  required 
for  the  tail  of  the  expansion  wave  to  reach  the  downstream  end  of  the  charge  tube  and 
the  value  of  x  is  the  effective  distance  traveled  by  the  wave  tail.  In  Becker's  evaluation 
of  the  boundary-layer  parameters,  he  assumed  that 

u/ui  =  (y/6)i/7  .23'i 


p  -Pw  _  u 
Pi  -  Pw  ui 


(24) 
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Cfj  =  0.045  Reg-^/'* 


(25) 


and 


Cfi/Cf  =  Fc 


(26) 


The  accuracy  of  this  choice  of  parameters  is  optimum  for  values  of  Reg  of  about  100,000 
but  deteriorates  as  the  Reynolds  number  increases. 


4.0  SKIN-FRICTION  COEFFICIENT 


One  widely  used  expression  for  incompressible  skin  friction  which  correlates  well 
with  experimental  data  over  a  wide  range  of  Reynolds  numbers  is  that  of  von  Kirm4n 
and  Schoenherr  (Ref.  6) 


(0.242)2 

“  {Jog  Regj  +  1.1696)  {log  Regj  +  0.3010) 


(27) 


which  can  be  integrated  to  give  the  familiar 

Cf.*/2  =  o.242//og(2Re9i)  (28) 


Another  expression  (Ref  7)  is  that  based  on  Coles'  law  of  the  wall  and  law  of  the  wake 
K  (2/Cfj)^/2  =  +  0.5  ^  (Q./2)  +  kC  +  2n  (29) 


where  the  constants  k  and  C  are  0.41  and  5.0,  respectively,  and  fl  is  a  function  of  Reynolds 
number  and  pressure  gradient.  If  the  laminar  sublayer  is  neglected  and  the  wake  function 
is  represented  by  a  sine2  distribution,  integration  of  the  profile  gives 


6i*  1  +  n 


(30) 


and 


T"T‘^ »i) 


Equations  (30)  and  (31)  must  be  used  with  Eq.  (29)  to  determine  Cfj  as  a  function  of 
Regj.  It  may  be  noted  that  for  an  earlier  version  of  the  wake  distribution  (Ref.  8),  the 
coefficients  of  fl  and  112  jn  gq.  (31)  3  2  and  1.522,  respectively. 
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In  order  to  use  Eq.  (29),  values  of  11  must  be  known  as  a  function  of  Reynolds 
number,  even  for  the  zero  pressure  gradient  condition  assumed  herein.  For  the  wake 
distribution  used  in  Ref.  8,  Coles  found  that  11  seemed  to  have  a  constant  value  of  0.55 
for  values  of  Re^j  greater  than  about  6,000.  For  the  sine^  distribution  used  to  obtain 
Eq.'  (31),  the  value  of  11  must  be  increased  slightly  to  about  0.56  at  Reo-  =  5,000  and 
further  to  about  0.58  at  Re^j  =  29,000  in  order  to  match  the  tabulated  values  of  Cfj 
therein.  Even  if  Eq.  (29)  could  be  put  into  a  form  which  could  be  integrated,  the  question 
arises  as  to  how  n  varies  with  Reynolds  number  to  values  about  three  orders  of  magnitude 
higher  than  that  considered  by  Coles.  Such  large  Reynolds  numbers  would  be  encountered 
in  a  large  Ludwieg  tube.  Even  the  high  Reynolds  numbers  of  the  experimental  values 
of  Kempf  (Ref.  9)  are  about  two  orders  of  magnitude  too  low  for  comparison  with  a 
large  Ludwieg  tube. 

If  it  is  assumed  that  both  Eqs.  (27)  and  (29)  give  identical  results  and  Eq.  (31) 
is  used  to  relate  Re^j  with  Re^,  values  of  11  can  be  calculated  by  an  iterative  method. 
The  results  are  shown  in  Fig.  3  and  indicate  an  increasing  value  of  11  with  increasing 
Reynolds  number. 


A  third  expression  for  skin-friction  coefficient  used  in  Ref.  10*  is 

_  _ 0.0773 _ 

(/og  Reej  + 4.561)  (tog  Rcflj- 0.546) 

*Constant  4.561  in  Eq.  (32)  was  4.563  in  Eq.  (70)  of  Ref.  10. 


(32) 
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which  can  be  integrated  to  give 

^  _ 0.0773 _ 

(to^  Re8j  + 4.0895)  (to^Reflj- 0.9431) 

Equations  (32)  and  (33)  have  been  found  to  correlate  well  with  available  data.  Equation 
(32)  agrees  almost  exactly  with  Coles'  tabulated  values  for  Reo-  from  4,000  to  29,000. 
It  gives  values  about  two  percent  less  than  Eq.  (27)  at  Re^j  =  500,  the  same  as  Eq. 
(27)  at  Re®  j  =  23,300,  and  about  six  percent  greater  at  Re^  j  =  20,000,000.  If  11  is  calculated 
from  Eqs.  (32),  (29),  and  (31),  the  other  curve  in  Fig.  3  is  obtained.  In  this  case,  11 
has  a  maximum  value  of  about  0.S88S  at  Reg  about  52,000  and  decreases  at  higher 
Reynolds  numbers.  Although  the  value  of  Cfj  from  Eq.  (33)  is  less  than  two  percent 
greater  than  that  from  Eq.  (28)  at  the  highest  Reynolds  number  of  Kempf's  data,  it  is 
believed  that  further  extrapolation  by  Eq.  (32)  is  better  than  that  by  Eq,  (27)  and  therefore 
Eqs.  (32)  and  (33)  are  used  hereinafter.  For  obtaining  the  ratio  of  dj/S,  the  combination 
of  Eqs.  (32),  (29),  and  (31)  is  used. 

In  Fig.  4,  the  limited  range  of  application  of  Eq.  (25)  is  clearly  shown  in  comparison 
with  Eq.  (32)  or  even  (27). 


Figure  4.  Comparison  of  incompressible  skin-friction  relations. 


Conversion  of  the  incompressible  skin  friction  coefficient  to  the  compressible  value 
with  heat  transfer  uses  the  concept  that  the  value  FcCf  (or  FcCp)  is  a  function  of  Fr^ 
Regg  in  the  same  manner  that  Cf.  (or  Cpj)  is  a  function  of  Regj.  The  factor  Fc  used 
herein  is  that  used  by  Spalding-Chi  (Ref.  11)  and  Van  Driest  (Ref.  12) 

Fc  =  (P/Pl)^^^  d(u/ui)]'^  (34) 


13 


AEDC-TR-75-118 


where  p/pj  =  Ti/T  within  the  boundary  layer  inasmuch  as  the  static  pressure  in  the 
boundary  layer  is  assumed  to  be  constant.  The  temperature  distribution  in  the  boundary 
layer  is  assumed  to  be  the  quadratic 

T  ^  Tfl,  (Taw  'Tw)  u  ^aw 

t7  ^  IT  Ti  ■  \  Ti ' 


where 


Taw/Ti  =  1  +  r(7  -  l)Mi2/2 

After  substitution  and  integration,  Eq.  (34)  becomes 

_  (Taw/Ti)-1 
(sin-1  a  +  sin-1 


where 


a 


Taw  Tw  ■  2T | 

[Taw  +  T„)2  -4Ti  Twl^/^ 


(36) 


(37) 


(38) 


and 


Taw  •  Tw 

,  ^  "  [(Taw+Tw)2-4T,  Tw]>/2  ’ 

The  factor  Fr^  used  herein  is  that  suggested  by  Van  Driest 

Frj  =  Pi/Pw  (40) 

and  Sutherland's  viscosity  law  is  used.  These  factors,  Fc  and  Frj,  have  been  used  for 
correlation  at  supersonic  speeds  with  good  results  and  should  be  satisfactory  at  the  relatively 
low  speed  and  heat  transfer  within  the  Ludwieg  tube. 


An  additional  factor  which  needs  to  be  established  is  the  velocity  distribution  in 
the  boundary  layer  so  that  Eqs.  (9),  (10),  and  (13)  can  be  evaluated.  For  simplicity, 
the  power-law  velocity  distribution  is  assumed 


u 

ui 


(41) 


14 


AEDC-TR-75-11B 


where  n  is  a  function  of  Reynolds  number  and  increases  with  increasing  Reynolds  number. 
If  a  kinematic  momentum  thickness  is  defined  as 


(42) 


then 


0k /6  =  n/(n2  +  3n  +  2)  (43) 

The  best  correlation  with  data  was  found  if  it  were  assumed  that  the  value  of  n  is  a  function 
of  Reynolds  number  based  on  the  actual  boundary  thickness,  not  corrected  by  Fr^,  and 
dk/S  given  by  Eq.  (31)  with  11  evaluated  from  Eqs.  (29)  and  (32)  with  dk  used  instead 
of  0i.  Then 


=  ^  Sa 

"  2  j0k 


3  + 


Values  of  n  evaluated  in  this  manner  are  shown  in  Fig.  S. 


(44) 


Figure  5.  Variation  of  n  with  Reynoids  number. 
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If  a  kinematic  displacement  thickness  is  defined  as 

('  -  TT-)  («) 

other  definitions  of  n  could  be  made 

n  =  -  1  (46) 

or 

”  =  iWJ  w’) 

The  three  definitions  of  n  give  the  same  value  at  an  Re^  of  about  50,000,  but  the  valiie 
defined  by  £q.  (47)  is  about  six  percent  higher  than  that  defined  by  Eq.  (44)  at  Res 
of  10^.  Inasmuch  as  11  is  evaluated  from  the  ratio  dk/fi>  it  seems  more  logical  to  also 
evaluate  n  from  the  same  ratio,  and,  in  fact,  Eq.  (44)  does  seem  to  produce  better 
correlation  with  experimental  data. 

5.0  CALCULATION  PROCEDURE 

The  computer  program  for  calculating  the  boundary-layer  growth  is  given  in  Appendix 
A.  A  design  Mach  number  in  the  Ludwieg  tube  is  assumed  and  used  to  obtain  an  inviscid 
sonic  area  from  Eq.  (1).  It  is  assumed  that  the  boundary  layer  at  the  sonic  area  location 
is  related  to  the  boundary  layer  at  the  end  of  the  Ludwieg  tube  by  streamlines  within 
the  boundary  layer.  Along  each  of  these  streamlines,  constant  total  pressure  and 
temperature  (enthalpy)  are  assumed  which  vary  from  streamline  to  streamline  as  determined 
from  the  assumed  distributions  of  velocity  and  temperature  within  the  boundary  layer 
at  the  end  of  the  Ludwieg  tube.  The  resulting  distributions  can  be  integrated  to  give 
a  displacement  thickness  which  reduces  the  effective  sonic  area  as  well  as  the  area  in 
the  Ludwieg  tube.  The  Mach  number  in  the  Ludwieg  tube  thereby  varies  slightly  with 
time  inasmuch  as  the  "viscid"  area  ratio  differs  from  that  assumed  initially. 

The  earliest  time  for  which  calculations  are  made  is  the  "start"  time  previously  defined 
in  Fig.  2.  As  a  first  approximation,  the  design  Mach  number  is  used  to  determine  the 
pressure,  temperature,  velocities,  and  corresponding  unit  Reynolds  number.  Because  Cp 
is  based  upon  Re^^  and  n  is  based  upon  Reg,  successive  approximations  are  made  until 
a  consistent  set  of  values  are  determined  to  solve  Eq.  (21)  for  5  and  to  obtain  the 
corresponding  values  of  displacement  thickness  both  at  the  end  of  the  Ludwieg  tube  and 
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at  the  sonic  area  location.  From  the  latter,  a  second  approximation  for  Mach  number 
is  obtained  and  the  process  is  repeated.  Usually  about  five  complete  approximations  are 
needed  to  make  the  Mach  number  consistent  with  the  "viscid”  area  ratio. 

From  the  final  calculation  for  the  "start"  time,  the  effective  location  of  the  origin 
of  the  expansion  wave  is  determined  for  use  in  the  calculations  for  later  times.  For  the 
later  times,  calculations  are  made  both  at  the  tail  of  the  wave  and  at  the  end  of  the 
charge  tube.  In  the  solution  of  the  momentum  equation,  it  was  assumed  that  certain 
variables  were  constant,  but,  in  actuality,  there  is  some  variation  along  the  tube.  Therefore, 
at  the  end  of  the  charge  tube,  the  denominator  of  Eq.  (20)  is  the  arithmetic  average 
between  that  at  the  tail  of  the  wave  and  that  at  the  end  of  the  charge  tube.  Approximations 
are  made  until  the  Mach  numbers  at  the  tail  of  the  wave  and  at  the  end  of  the  tube 
are  consistent  with  the  "viscid"  area  ratios  at  the  corresponding  locations.  For  each 
approximation,  it  is  assumed  that  the  Mach  number  at  the  tail  of  the  wave  determines 
the  stagnation  pressure  in  the  tube,  but  the  Mach  number  at  the  end  of  the  tube  determines 
the  static  pressure  at  the  end  of  the  tube.  Again,  about  five  or  six  approximations  are 
needed  to  achieve  consistency. 

Obviously,  the  maximum  value  of  5  is  the  radius  of  the  charge  tube.  The  calculated 
time  at  which  this  condition  first  occurs  is  the  maximum  time  for  the  particular  Ludwieg 
tube  design  even  though  the  reflected  expansion  wave  may  not  have  returned  to  the 
downstream  end  of  the  tube.  If  there  were  no  boundary  layer,  the  pressures  within  the 
tube  would  be  constant  during  each  excursion  of  the  expansion  wave  up  and  back  down 
the  tube.  In  actuality,  the  static  pressure  at  the  downstream  end  of  the  tube  decreases 
with  time.  This  decrease  becomes  greater  as  the  design  Mach  number  in  the  tube  is 
increased.  A  phenomenon,  reported  by  Piltz  in  Ref.  13  and  qualitatively  supported  by 
calculations  made  by  Piltz  and  by  the  method  described  herein,  is  that  the  stagnation 
pressure  which  initially  decreases  slightly  with  time  may  increase  slightly  at  later  times 
at  design  Mach  numbers  above  about  0.2  but  continues  to  decrease  at  lower  design  Mach 
numbers.  The  magnitude  of  these  decreases  and  increases  is  a  fraction  of  one  percent 
of  the  theoretical  inviscid  values.  After  the  boundary-layer  thickness  becomes  equal  to 
the  tube  radius,  the  stagnation  pressure  decreases  more  rapidly  at  first  and  less  rapidly 
later  as  the  boundary  layer  adjusts  to  the  velocity  profile  of  fully  developed  tube  flow. 
Such  behavior  can  be  seen  only  if  the  Ludwieg  tube  is  sufficiently  long.  From  a  practical 
standpoint,  it  would  not  be  economical  to  build  a  Ludwieg  tube  with  a  run  time  greater 
than  that  which  would  allow  the  boundary-layer  thickness  to  become  nearly  equal  to 
the  tube  radius. 
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6.0  COMPARISON  WITH  EXPERIMENT 

The  comparison  shown  in  Fig.  1,  which  was  taken  from  Ref.  4,  was  made  before 
the  method  described  herein  was  developed  to  its  present  form,  primarily  in  the  realm 
of  the  variation  of  n  with  Reynolds  number  and  the  incorporation  of  the  stream  tube 
method  of  changing  the  effective  sonic  area  with  time.  A  comparison  is  made  in  Fig. 
6  of  the  same  experimental  values  with  values  calculated  as  described  herein.  Two  sets 
of  calculated  values  are  shown  to  illustrate  the  sensitivity  of  the  method  to  changes  in 
start  time  and  run  time.  One  set  of  values  uses  the  same  start  time  of  0.032  sec  for 
each  pressure  level  together  with  run  times  of  0.060  and  0.110  sec.  Curves  drawn  through 
these  values  indicate  a  greater  influence  of  Reynolds  number  on  boundary-layer  thickness 
than  that  shown  in  Fig.  1  and,  therefore,  agree  somewhat  better  with  the  experimental 
values.  Moreover,  the  bands  of  the  experimental  data  indicate  the  inaccuracies  and  spread 
of  the  data  from  many  runs.  The  times  for  the  second  set  of  calculated  points  in  Fig. 
6  were  selected  so  that  the  points  for  the  lower  and  higher  charge  pressures  lie  more 
nearly  in  the  center  of  the  experimental  bands  in  the  manner  of  the  points  for  the  medium 
pressure.  Only  small  changes  in  the  times  were  necessary  to  produce  the  changes  in 
boundary-layer  thichness.  These  slight  time  variations  which  produce  improved  agreement 
are  well  within  the  experimental  uncertainty. 


Range  of 


Figure  6.  Comparison  of  calculated  and  experimental  results  for 
11.75-in.-diam  Ludwieg  tube,  revision  of  Figure  1. 
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The  discussion  above  illustrates  what  the  dangers  may  be  in  attempting  to  extrapolate 
experimental  data  to  other  sizes  and  operating  conditions  without  an  adequate  theoretical 
method  as  a  backup.  The  length  of  the  Ludwieg  tube  used  was  insufficient  to  give  values 
of  6  greater  than  about  one-half  of  its  radius.  The  agreement  between  calculated  and 
experimental  results  in  the  7.6-percent  scale  pilot  tunnel  indicate  that  the  present  method 
would  be  adequate  for  application  to  the  design  of  a  full-scale  facility.  Data  from  charge 
tubes  of  11.75-  and  13.94-in.  diameter  were  obtained. 

Relatively  few  boundary-layer  data  were  obtained  at  the  downstream  end  of  the  larger 
charge  tube.  These  are  compared  with  calculated  values  in  Figs.  7  and  8.  The  full 
boundary-layer  thickness,  6,  (Fig.  7)  is  predicted  quite  accurately.  Close  to  the  wall, 
however,  the  mass-flow  profile  (Fig.  8)  is  underestimated  by  the  calculations,  but  beyond 
about  0.5  in.  the  agreement  becomes  much  better.  Because  of  the  profile  deviations,  the 
values  of  displacement  or  momentum  thickness  are  not  quite  as  accurately  predicted  as 
the  full  thickness. 


Figure  7.  Comparison  of  calculated  and  experimental 
results  for  13.94-in.-diam  Ludwieg  tube. 
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Figure  8.  Calculated  and  experimental  mass  flux  profiles 
in  the  13.94-in.-diam  Ludwieg  tube. 

For  the  range  of  Mach  numbers  under  consideration,  an  error  of  O.S  percent  in 
measured  pressure  ratio  will  result  in  an  error  of  nearly  five  percent  in  Mach  number 
near  the  edge  of  the  boundary  layer  and  an  error  of  nearly  ten  percent  in  Mach  number 
at  a  point  in  the  boundary  layer  where  the  Mach  number  is  about  70  percent  of  the 
free-stream  value.  Thus,  the  deviations  between  calculated  and  experimental  profiles  are 
within  the  accuracy  of  the  experimental  values. 

7.0  CONCLUDING  REMARKS 

A  method  is  presented  for  calculating  the  time-dependent  growth  of  the  boundary 
layer  at  the  downstream  end  of  a  Ludwieg  tube.  The  method  consists  mainly  of  several 
modifications  to  Becker's  method  which  had  been  found  to  be  inadequate  for  use  at  high 
Reynolds  numbers.  Calculations  made  by  the  modified  and  improved  method  have  been 
found  to  agree  quite  satisfactorily  with  experimental  data  obtained  in  a  Ludwieg  tube 
used  to  drive  a  small  transonic  tunnel.  Utilizing  the  program  developed  in  this  study, 
very  reliable  predictions  of  the  boundary-layer  growth  in  the  charge  tube  for  wide  range 
performance  are  possible. 
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APPENDIX  A 
COMPUTER  PROGRAM 

This  Appendix  contains  a  program  listing  together  with  a  sample  output  of  calculated 
results.  The  computer  program  is  written  in  double  precision  Fortran  IV  for  use  with 
the  IBM  370/165  Computer.  It  closely  follows  the  description  given  in  Calculation 
Procedure.  Numerical  integration  is  accomplished  through  the  use  of  a  1 6-point  Gaussian 
formula  for  the  interval  0  to  1.  To  avoid  the  problem  of  infinite  slopes  when  y/6  is 
the  independent  variable,  the  velocity  ratio  is  made  the  independent  variable,  because, 
from  Eq.  (41) 

y/6  =  (u/ui)"  (A-1) 


and 


d(y/6)  =  n(y/8)d(u/ui  )/(u/ui )  (A-2) 

In  the  program  u/ui ,  written  in  the  output  as  U/UE,  is  the  variable  Z(K)  and  the  D(K)'s 
are  the  corresponding  weighting  factors  for  the  Gaussian  integration.  One  subroutine, 
FMR,  is  used  to  determine  the  Mach  number  for  a  radius  ratio. 

Four  cards  supply  the  input  data  for  a  particular  problem.  The  first  card  contains 
the  title  (ITLE)  information  in  columns  2  through  12.  On  the  other  three  cards,  the 
format  allows  ten  columns  for  each  variable. 

Second  Card 


Input 

Columns 

GAM 

1  -  10 

Ratio  of  specific  ^ats,  y 

AR 

11-20 

Gas  constant,  ft^/sec^R 

ZO 

21  -  30 

Compressibility  factor,  1 

RO 

31-40 

Recovery  factor,  r 

Vise 

41  -  50 

Constant  in  viscosity  law 

VISM 

51  -  60 

Constant  in  viscosity  law,  viscosity  = 

VISC(T)i-5/(T  +  VISM) 
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Third  Card 


AKAT 

1  -  10 

Variable  used  if  more  than  one  problem  is  input 
at  same  time,  see  card  A-290 

CMACH 

11-20 

Design  Mach  number,  if  not  specified,  maximum 
value  is  calculated  from  the  radius  ratio, 

PP/RSTAR 

PP 

21  -  30 

Radius  (in.)  of  Ludwieg  tube 

RSTAR 

31  -  40 

Effective  radius  (in.)  of  test  section  or,  if 

CMACH  =  0,  of  sonic  area;  if  RSTAR  =  0,  RSTAR 
set  equal  to  sonic  radius  calculated  from  CMACH 

PPQ 

41-50 

Charge  pressure,  P4,  psia 

TQ 

51  -  60 

Charge  temperature,  Rankine 

TMIN 

62  -  70 

Start  time,  sec 

BLD 

71  -  80 

If  zero,  boundary-layer  profiles  are  not  printed 

Fourth  Card 

TM(K) 

10  for  each 
value  of  K 

Times,  seconds,  at  which  calculations  are 
desired;  maximum  value  of  K  is  7;  problem  is 
terminated  when  TM(K)  =  0 

Output  values,  if  not  otherwise  obvious 

RE/IN  Reynolds  number  per  inch  in  freestream 

RTHI  Incompressible  Reynolds  number  based  on  momentum  thickness 

FRD  Rfl./R,^ 

FC  Cfi/Cf 

KCFI  1000  Cfj 

KCDI  1000  Cpi 

KCD  1000  Cf 

TH  di ,  in. 
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Y 
M 

D*2-D 

DELTA* 

DEL*A 

H 

R 

PO 

POO 

PE 

PEO 

TO 

TAU*V 

XSTAR 

HRH 

TAO 

DELAY 

TIME 

V 

HSUM/VO 

HSUM/VE 

THROAT 

TUBE 

NET  TIME 


R-DEL*A 

Mach  number  at  end  of  Ludwieg  tube 

f® 

I  (1  -  pu/piui)dy 
0 

6i,  in. 

6*,  in. 

6i/0i 

Tube  radius,  r,  in. 

Stagnation  pressure,  psia,  at  net  time 
Stagnation  pressure,  inviscid,  at  time  0 
Static  pressure,  psia,  in  freestream 
Static  pressure,  inviscid,  at  time  0 
Stagnation  temperature  at  net  time 

Distance  to  concentrated  expansion  wave  from  end  of  tube 

Distance  to  end  of  tube  from  effective  origin  of  expansion 
wave 

Distance  to  head  of  expansion  wave  from  end  of  tube 

Time  required  for  head  of  expansion  wave  to  reach  end  of  tube 
from  effective  origin 

Velocity  of  concentrated  wave,  fps 

(p*  =  Si)/0iUi  at  tail  of  wave 

At  end  of  tube 

Values  at  RSTAR  calculated  from  streamtube  from  end  of 
Ludwieg  tube  to  sonic  area 

Length  of  Ludwieg  tube  for  NET  TIME  for  expansion  wave  travel 
from  downstream  end  of  tube  to  upstream  end  and  return 

Same  as  input  TM(K);  useful  run  time.  Fig.  2,  is  NET  TIME 
minus  START  TIME 
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C  PRU6RAM  LUDt>UG 

C  dOUNOANY  LAYER  GHOdTn  IN  lOOdIEG  ludb 

Implicit  real»8  lA-rno-^i 

Common  /GG/  G1 »G3tGG«G6tGAtRGA 

OIMENblON  Z(l«)t  D(lb)t  ITlE(3)>  1M(SI 

OnEmI.U'U 

ZCRObO.OO*U 

0*l)s. 01357622910*0 

U(2)b. 03112616200*0 

0<3)«. 04^57925580*0 

0|4)>. 06231448560*0 

0(5)>. 07479799440*0 

0I6)>. 08457825970*0 

0(714.09130170750*0 

0(8)4,09472530520*0 

2(1)4.00329953250*0 

2(2)4.02771248850*0 

2(3)4.06718439880*0 

2(4)4.12229779580*0 

2(5)4.19106187760*0 

2 (6) 4.27099161120*0 

2(7)4.35919622460*0 

2(8)4.45249374510*0 

00  1  jB9tl6 

d(J)4U(17-J) 

1  2(J)4l.O*0-2(17>J) 

2  HEAU  {5,30,E90b29)  I TLL 

3  REAO  (5*31 f EN0429)  GAM.ANtZO.ROtOISCtVISM 

C  FOR  GAMMA4l,4,  61*2. 5,  6341,5.  G54i/6*  6645/6,  6741, 2*  6840.2 

C  GA43,  RGA4123*  GM40.4,  GP43.5 

6M4GaM>1.0*0 
6141,0*0/69 
684.50*0469 
6741,0*0*66 
6641,0*0/67 
65468*66 
GA4G)4G7 
63>.50*04GA 
6P46AM461 
RGA46M46b 

4  READ  (5,31,END429)  A«AT .CMACM.PP.RSrAR.PPO. lO.TMIN.dLO 

C  RSTAR(HlRT) 460.5551605,  RSTAR(PlLOT  MlRT) 44.62193583 

READ  (5*31,L90429)  ( TMIK) .KbI.R) 

MC4l,D*2*C9ACH 

IRST4l,0*14RSTAR 

IF  ((MC.£Q.O).ANO.(IRST.Ne.O))  CMACHbEmR  (PP/M5TAR  ,-1^) 
BMMbCMACH 

OMACMbCMACR  _ 

■RITE  (6,40)  ITLE.PPU.TO.THIN.CMACM 
■RITE  (6,33) 

YSTARbPP^OSORT (CMACM)/(G6*G5*CMACH442)44G3 
IF  (9C.EO.0)  Y5TAR4R5TAR 
IF  (IRST.EO.O)  RSTAH4YSTAR 
XPT4RSTAR/YSTAR  _  _ 


i 

i 

3 

6 

7 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

21 

22 

23 

24 

25 

26 

27 

28 

29 

iU 

31 

32 

33 

35 

36 

3f 

38 

39 

40 

41 

43 

44 

43 

46 

47 

48 

40 

SO 

5l 

52 

53 

S4 

26 


AEDC-TR-75-118 


S)MACM»FKf<(X«Ti-l) 

4 

SB 

*0>USaHT(G*M>AR*ro) 

A 

5b 

Tx»IO 

A 

“57 

>IMU>VlSC*TW*DSURT(tM)/tTN*VISH) 

A 

58 

IQaO 

a“ 

59 

TIMiTMIN 

A 

6U 

CAPI-t&SU*# 

"a 

61 

TAO>12.O*0*TIM*Ad 

A 

62 

[>LB«0>OU«0 

A 

63 

1>LC«0.0I)*0 

A 

6* 

•«*sO 

r 

65 

1PP»0 

A 

66 

HAaHA»l 

A 

6?' 

lfiS»l.U*U/ ( 1 .0»0»6B»OHACH) 

A 

68 

OatTmbb»OMACH»»2 

A 

69 

STH>l.D*0/(l.D*a*OdEr) 

A 

70 

Ih»TRS/STH 

A 

71 

T0«TR»TQ 

A 

72 

pps«ppu»tm»»sp 

A 

73 

IF  ( (tU.tU.O) .and. (HA.EU.1I )  PS0«PP5 

A 

74 

Hh0>1*a.U«V*PPS/20/AH/10 

A 

75 

lE«TO»brH 

A 

76 

£MU»VISC»TE«3SQHT (TE)X(TE»VISHI 

A' 

“7> 

tANaTE*(l.D*0*H0*0SEI) 

A 

78 

liHOEmRHO^STP^^Gl 

A 

"79 

AEeOSQHT (gahoar*te ) 

A 

80 

ve«omach«ae 

A 

ST 

vo>omach*ae 

A 

82 

AT»AE“V0 

A 

“6J 

IF  (TIH.GT.TH1N)  GO  TO  7 

A 

84 

OSEC«TM1n»AT/(AQ-AT) 

A 

~S5 

*ST-AQ*UbEC 

A 

86 

AV«(AO*2.0«0>AT)/i.U*0 

A 

8f 

TSEC»TIM»OSEC 

A 

8B 

TAUV>12.U*U*TS£C*(A)/>A() 

A 

89 

«EO»RHOt»VO/EMU/l2.G*0 

A 

90 

HaT'O 

a" 

91 

GO  TO  V 

A 

92 

IPP»0 

A 

93 

TaUV>12.g*0*( TSEC*Ay-XST) 

A 

94 

PPt*PPS/  ( 1 .0*  0«GB««MN»«<:)  *»GP 

A 

95 

IF  (  (lO.EU.O)  .and.  (HA.EU.D)  PEU«PPE 

A 

96 

BET«GB»hhN«*2 

A 

97 

STR>1.0*0/(1.0*0*Ptr) 

A 

98 

TE«T0»BTH 

A 

99 

Iah«TE»(1.D»0»HO»bET) 

A 

100 

HHOt«HMO*STH»*Gl 

A 

id' 

eM0»Vl»C»TE»D5QRT (1E)/(Tt»VISH) 

A 

102 

VE«NHN*(}SURT  (GAM*AR*(E) 

A"" 

ioJ 

KEO«RHO£»VE/EHU/12.0»0 

A 

104 

i>G>T.M/IE 

A 

iOB 

0a«Ta«-Tm 

A 

106 

i»B»TAti-TE 

A 

i07 

OK.Da/TE 

A 

108 
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UNsOy/Tb  *  109 

DC>USU**T  1IIA*3A**.03*U«1«>U3)  A  UU 

OF^QAN&IM  (Ud*T9-Tbl/UCI  A  Ill 

[>b>UA*<^l9(DA/ULI  *  112 

ii»ua/(ur*uf.)/(DF<iit)  A  Hi 

FrtUaCHu/MMJ  A  llA 

XSTAI««TAUV/<1.U*U*1.JO*0*AV/vE)  a  lib 

FC»TP/It  A  116 

9XI>raO*rib3*XSTA>I/FC  A  117 

IF  (RXI ,LT.1.U*2I  so  to  26  A  116 

ABAUL001u(9Xll-l.a3*U  A  119 

>4THI«U.U*AU*0*MXI/Ad*A2  A 

MTlIiHlHl  A  Ul 

M0bL>1.0*l*)lTHI  A  122 

10  IF  (RTHl.LT. 1.0*1)  00  TO  27  A  120 

9THG«0L0(ilU14THl>  A  12A 

RTtbBOuOblU (MTll )  A  126 

I.DlBO.UiHb6D*0/(i<T>lO*A.U9950*0>/(yrHb-0.94Jlu*UI  A  126 

CFI>O.UJb660*0/t>«T10*A.6610*0)/t>lT10-U.!>46|i*U)  A  127 

XCF«.a10«0/OSORT (CFI)  A  i2a 

11  (;2Bl.0*0*CAPl  A  129 

CJb2.0*0*CAPI<(3. 176979720*0*1. 60*0*CAPII  A  139 

Cl«C2-L3/xCF  A  131 

FACFaXCF*OLOGICll -1.166«016810*u-2.O*U»CAPI'2. JU26bbU9JO*0*HTIb  A  132 

FPCP«(XCF-3. 178979720*9-3. 0*0*CAPl)/XCF/Cl-2. 0*9  A  13J 

CAPI»CAP|-FXCF/FPCP  A  13A 

IF  (DA6SIFXCF) .b1 .1.0-8)  bO  TO  11  A  13b 

U0TIbXCF/C1  a  136 

XNB.SO*d*(UOTl*0SU'lT(noTI*(0UTl-6.0*0)*l.O*U)-3.0*U)  A  137 

CObC01*1E/TP  a  138 

bUHAB0.0O*0  A  139 

SUH8b0.00*0  a  140 

bUt1L«0.UO*0  A  I4l 

bUHOBO.Oj*0  A  142 

bOHEBU.0j*0  A  143 

bMFB0.U0*0  A  144 

bHGBQ.91)*0  A  14b 

bMHB0.90*0  A  146 

1)0  12  XBlfl6  A  147 

JNaZ(K)**A9  A  l48 

IObO«*2IX)*(3K-2(K|40N)  a  i49 

OObO(K)*AN409  a  159 

AnUBOO/l-<  A  Ibl 

dUOBADu'ilO  A  152 

C0DbaI>O*9N  a  153 

>}DO«8D04ol*  A  154 

5UMABbu»lA*AU0  A  l3b 

bUMbBS061*8DJ  A  l3b 

bUMLaSurtC*COJ  A  137 

L0UBAUO/2IX)  A  Ibd 

•iOO*EOIJ*JN  A  159 

5UHb>5uM) *EUO  A  160 

5MhbSMH*603  a  161 

12  50MOb50Mj*OOJ  a  1o2 
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0OT»l.U«0/(SUMA-5U'<tt» 

A  163 

US00>1.1>*0-5JMA 

A  166 

OSM«0«50*0-5JMC 

A  165 

THHbSUHC-5UH9 

A  166 

801}«SUME-1.D*0 

A  167 

HMO«,5U*0-SNH 

A  16b 

If  (IPP.aT.U)  GO  TO  13 

A  169 

M»OSOO»OOT 

A  170 

OOTHsUOT 

A  171 

MRH«900«JOT 

A  172 

13 

IF  (KAT.EQ.0I  30  TO  U 

A  173 

HOV£«(H«hKH)/VE 

A  174 

ASTAR«TAJV/(1.0*U*.50*0*AV*(HOVO*HOVE) ) 

A  175 

60  To  lb 

A  176 

1* 

H0V0«(H^HHH)/V0 

A  177 

XSTAR>TAuV/(1.0*U*AV*HOVO} 

A  17B 

15 

1H«XSTAR*C0 

A  179 

()OR«OOTR*TH/PP 

A  180 

IF  (DOR. LT. 1.0*0)  SO  TO  16 

A  TbI 

DOR>1.U*0 

A  182 

Tm»PP/UOTR 

A  183 

16 

ushoo*dsjd-oor*osm 

A  184 

ooth«i.o*o/(i.o*o/oot-tm>i«oor) 

A  18b 

HRbOSROO*OOTR 

A  186 

riRH«<WUl>*O0R*RHD)*3UrR 

A  187 

IF  (D*b5(H-HR) .LT.5.0-7)  60  TO  17 

A  188 

HaHR 

A  189 

ipp»j*ipp 

A  190 

GO  TO  13 

A  191 

17 

OElTa*OOTR«TR 

A  192 

IF  (DOR. 60.1.0*0)  OtUTA'PP 

A  193 

«OLT«RtO*OELTA 

A  194 

HTHX.FR0*RDLT/'D0T 

A  195 

HTIl.RDLT/OOTI 

A  196 

0ELST»MR«TH 

A  197 

IF  (DABS(1.0*0-RThX/MTm1 ) .UT. 1.0-6)  GO  TO  Id 

A  198 

HTHIsRTMX 

'A  1^ 

ROELbROLT 

A  200 

GO  TO  10 

A  201 

18 

IF  (KAT.GT.O)  go  TO  19 

A  202 

OybOSORT (PP*(PP-d.O*OBOtLbT)  ) 

A  263 

KaT*! 

A  204 

HObH 

a‘ToT 

HRu«HR^ 

A  206 

XNO»XN 

A  207 

GO  TO  6 

A  208 

19 

COLSTB<:.0*0*OELST/'(OSORt  ( 1 .0*0-2.0*0*oELST/PP)  *  1  .0*0 ) 

A  209' 

0lFFBDAaS(C0LST-(<'.U*0*OtU*0LC)/J.0*0) 

A  210 

IF  (D1FF.LT. 1.0-6)  GO  TO  21 

A  211 

FSbOSORT (PP*(PP-2.0*UBUtL6T) ) 

A  212 

TTRBBtT*SIR/Gb 

A  213 

OLCbDLO 

A  214 

OO  20  KBj.ib 

A  21§ 

ONbZ(K)»bX» 

A  216 

29 
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■■  TH»ow*i;rK)»(i)K‘-i(K)»o»iiT  Cm 

00>0(K)*XN*UM  X 

UBR«DSUaT(tH-11H*(TH-ZIK)**2)  )  A  <il9 

FDD-DO/UhR  a  220 

SmF«5MF»F00  *  2Z1 

20  bM(»»5M&»f00»l^  _ _ _ _  A  2?2 

FF6-DSI-MT  1(1  .0*0-O0H)  ••2*2.O»0»OOH*  (SMF^liORiSMS)  I  A  223 

0LB»CW.Sr  A  22* 

RX.pp*FFS/FbTAN  A  22» 

1)H«0Y»Ra/FS  .  .  _  A  226 

«MN«FHH<RX<-1I  a  22r 

0maCh»Fmr(0h»-1)  _6. 

bo  ID  b  A  22» 

21  KTAO>Reo«TAU  .  A  230 

OOL<D2LTA/TAO  «  231 

0SUL>CUL<>r/TAO  _  A  232 

TRhXiXRI  a  233 

ImaCH»S<XACH  _  _  _  A  234 

IF  ( (INST. EU. (I)  .UR.  INC. 23.01  I  FGSaFFIi  A23B 

IF  (BLO.cQ. 0.00*0)  GO  TO  23  A  236 

1FIM0D(IJ(2).'«E.U>  MHirt  <6««1)  4  23T 

•RITE  (6.34)  A  232 

•RITE  (6«39)  ZERUt2tROt2EROtD«>2ERa  A  23« 

00  22  ^■T^16  _  A  240 

UN»Z(K)**Xi'  '  A  241 

TtiL>UN*t)t:LTA  A  242 

TR«02*2U)»(DA-2(K)*UN)  a  243 

RH0Ub2(K)/Tn  a  244 

22  RRITE  (6.33)  UN.TBLtZtK) .TR.RHOU  A  249 

•RITE  16.39)  OnE.OELIa.O^.O^EiONE  A  249 

23  IF  ( (IR9I  .EO.O)  .OH.  InC.tO.U)  )  GO  Tu~29'  A'  24  7 

9mFb0.<<O*D  a  240 

9mg«0.UO*0  a  249 

THtTBG«)4rMACH4*2  A  290 

ITT»8tT43TR»H.D»0*IiiET)/rHET  A  291 

UO  24  Aal.16  A  292 

UN»2(K)»4XN  a  293 

1R»04*214)»<OX-21X)*()N)-  A  2S4 

UObO(K)4XN4U9  a  299 

UBK«D50KT(Th-TIT*irK-2(R)**2) )  A  296 

FODBDO/UtlR  A  297 

9MF«SM7*fOD  A  290 

24  3Mb»9Hb*fDD»JN  '  ”  A  299 

FGG«0SbRI ( (l.U*U-O3R)442*2.0*04UUR4<9nF-00R*9H6) )  A  260 

Thh«xRT4FF6/FG6  a  261 

TMAtH«FR«(TwR.-l)  A  262 

IF  (DAri9(TRR-lRRX) .LI .9.3-9)  GO  10  29  A  263 

IRRXbTHR  a  264 

GO  To  23  ■"  A  269 

25  OSBORbI. 0*0-11. 0*0-COLST/RP)/FGG  A  266 

UROhBl,i)*0-U.0*U-0OR)/FSG  A  267 

U2)>1.9U*0*I1.0*O-TAUV2Ta0)  A  26B 

ThL».50*0*TRS*|I9F.C»IAE*YO)*X9I)-X9T  A  260 

C111Kb20Uu.O*U*CDI  a  270 
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CDK>2000.U*U*C0  *  271 

CFIK>200U,D*0*CFI  *  272 

I1S2>DS00*1)ELTA  *  273 

PSR«PP5/t»S0  A  27A 

PER»PPE/PeO  A  27b 

•  «ITE  (b.*2)  REOtRTMA,FRD»FC»CFIKtCOIK.COI««TH  A  276 

FRITE  (6i36I  HA«FSiMMN«DELTA«OS2>DELbTtCOLStiHtXNtPP(PPSfPPE»TO«TA  A  277 

lOVtXSTARiHRH.OMACHtPbRtPERfTAOfQSECtAV  A  278 

FRITE  (6i*3)  XNO*HO«r1RO«HOVO>HOVEiDWT  A  279 

FRITE  (6|32I  SMACHtTRACHtDSB0R«U80RtTbL  A  200 

FRITE  (6i39)  RTAOtOOLiDbOLiTlM  A  201 

10»IQ*1  A  202 

TIM«T«11J)  a  203 

IF  (Tm.LE.THIN)  60  TO  26  A  204 

IF  ( (8UD.Ea.0.0D»0) .AND. (IQ.NE.b) )  GO  TO  b  A  205 

1F((BL0.NE. 0.00*0). AND. (HOO(ia<2).NE.O))  60  FOb  A  206 

FRITE  (6i40)  ITLEtPPRiTS.THINtCMACH  A  207 

FRITE  (6.33)  A  208 

60  To  b  A  209 

26  IF  (aKAT)  3.2.4  A  290 

27  FRITE  (6.37)  RTHI.ROEL  A  291 

GO  TO  29  A  292 

20  FRITE  (6.38)  RXI.FRU.REO.aSTaR  A  293 

29  STOP  A  294 

C  A  295 

30  format  (JA4)  A  296 

31  FORHaT  (8E10.0)  A  297 

J2  Format  (1h  .iox.< inviscio  throat  h>'.f0.6.*.  streamiube  throat  m>  a  296 

l'.F0.6.>.  0EL*A/RAU.«*.F8.&.'.  UELTA/RaO.*' .F8.6.  • .  TUBEH.  A  299 

2F7.2.»  FT.LG*  /  )  A  300 

33  format  (IH  .bX. 'QUADRATIC  TEMPERATURE  OlSTRlbUTION  bPAULOING-C  A  301 
IHI  reference  TEMPERAIUHE  VAN  ORIEST  REFERENCE  REYNOLDS  NO.*  /)  A  302 

34  format  (lH0.iax.<Y/OtLTA'.10X.'Y(INCH)'«llXV'U/UE'.13X.'T/TE'*llX.  A  303 

I'MASS  FLUX'  /)  A  304 

3b  Format  (ih  ibFiT.b  )  a  30b 

36  format  (IH  .'APPR'.U.*.  Y»*.F9.b.'.  H«*.F0.6.*.  DELTA-' ,F9,5.  A  306 

1'.  D*2-D>'.FB.b.'.  OELTAX'.FG.b.'.  OtLoA-' .FB.bt ' .  H>'.F9.6.  A  307 

2'.  N«'.F8.b  ./llX.'R-'.F9.bt'.  P0>'tF8.3.'.  PC>'tFa.3.'.  TO-'.  A  308 

3F0.3.'.  TAU-V>'.F9.2.'  IN..  XbTAR-' .F9.3. '  IN..  HRH-' .F9.6/11X.  A  309 

4*MAVE  LNu  M-'.F0.b.'.  PO/POO-' .F8.b. ' .  PE/PbU-' .Fb.b. ' .  TaO-'.  A  310 

bF9.2.'  IN..  DELAY  T IHE- ' .F9.6. '  SECONDS.  V-'.F9.3  )  A  311 

37  format  (IMl.'RTHI-'.Elb.T.'RDEL-'.Elb.T  )  A  312 

38  format  (lHli*RXl*'.Elb.7.'FH0-'.Elb.7.'RE0*'.Elb.7.*XSTAR-'.£15.7)  A  313 

39  Format  (ih  .lox. 'Reynolds  no.  (Tau)-' .ipei2.s.'.  DELrA/TAO-'.Ei2.b  a  3m 

1".  DEL*A/TA0>'.E12.b.*.  NET  T IME- '  .0PF9.6. '  SECONDS'  /  )  A  31b 

40  Format  (ihi.3A4.'  charge  pressure-'. f7. 2.'  psia.  charge  temperatu  a  316 

IRE-' .F7..:. '  R.  start  T1ME-'.F9.6.*  SECONDS.  DESIGN  M-'.F8.6  /)  A  317 

41  format  (IHOI  A  318 

42  format  (IHO.lOX.'RE/lN-'t  F9.0.'.  RTHI - * .F9. 0 . ' •  FRO-' .Fa. 6. ' .  A  319 

IFC-' .F8.b. I  I  aCFI-' .Fb.bi ' .  KCDI-< .F6.b.'.  KCO-* .F0.S. * .  TH-'.  A  320 

2F8.5  )  A  321 

43  format  (Ih  .lOX.'HAVL  EN3  N>'.F6.5.'.  Hc'.F9.6.'.  HRH- ' .F9.6. ' .  A  322 

1  HbUM2V0»'.F11.9.'.  MSUH2VE-'.F11.9.'.  HAVE  LGTH/TAO- * »F7.b  /)  A  323 

ENO  a  324- 
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Function  fin  (Rt**is)  ‘t  i 

C  10  obtain  MACH  NUMBtH  FN3M  RAUIU&  HAIlOi  IS«-1  IF  •*<!#  *1  IF  M>1  d  i 

Implicit  HeAL»b  (A-Mto-/i  a  3 

common  /Ub/  61 iGJtObtGbtOAtRGA  d  A 

IF  IRH.Gt. 1.0*0*6.0-101  60  TO  I  a  6 

fmr»i,o*ii  a  b 

IF  (RR.LT. I. 0*0-6. 3-10)  60  TO  b  d  7 

RETURN  a  6 

1  AH»(RH»»i)»*RCA  3  6 

IF  (Ib»A6.LT.-l.*!>0*W|  GO  TO  3  a  lu 

FHH>(l.U*U*Ib*DbURT((AR-l.D*b)/Gl)  lOGA  d  II 

60  To  3  d  12 

2  FMN»<bb»*G3/RH)»»2  d  13 

3  00  A  J«ltbO  3  lA 

CM>FmR**R6A  d  lb 

FH*CM»AB-6b-6b*FMH»*2  a  lb 

FP»H8A»(LM»Ad/FMH-FMH)  d  1/ 

FMH»FMH-FM/Fa  3  l8 

IF  IDAUSIFM) .LT. 1.0-9)  faO  TO  b  a  19 

A  CONTINUt  a  2U 

6  RSa (G6*b3*FMM**2) ••bJ/bbORT (FMN)  d  21 

MRITE  (bW)  RHiFMHtRb  d  22 

return  a  2J 

a  2A 

format  <IH  t'RAOIUS  HAT13»*.F12.10*»*  M» « .F 12. 1 9 .Tm»  M/R*»tF12. lUl  3  2b 

tNO  3  2o- 


gample  Input 

CARO  1 
ITLE 


PILOT  HINT 

CARO  2 

6AM 

AR 

20 

HO 

Vise 

VISM 

.1l.a 

CARD  3 

171b. 3b 

1. 

JL*«94 

2.2b97E-a 

._1^8,T2 

AKAT 

CMaCH 

PP 

RbTAH 

PPU 

TO 

THIN 

HLD 

1. 

CARO  A 

.26b 

6.96B7b 

A.b2I93be3 

102. 

b2U. 

.0A2 

TMtl) 

.ors 

TH(2) 

.lib 

TMUI 

.Ibb 

TMCAI 

TM(b( 

TM(b) 

IM(7) 
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PILOT  hirt  Charge  pressures  loz.oo  psut  chargc  temperature-  szo.oo  r>  start  time*  o.oazooo  seconos#  design  msq.zgsooo 
quadratic  temperature  distribution  SPAULOING'CHI  reference  TEMPEHaIURE  van  driest  reference  REYNOLDS  NO. 


y/oelta 

Y(INCH) 

U/UE 

T/TE 

MASS  FLUX 

0.0 

0.0 

0.0 

1.11084 

0.0 

0. 00082 

0.00038 

0.35920 

1.07403 

0.33444 

0.00406 

0.00191 

0.45249 

1.06392 

0.42531 

0.01526 

0.00717 

0.54751 

1.05339 

0.51976 

0.04551 

0.02139 

0.64080 

1.04282 

0.61449 

0.11142 

0.05236 

0.72901 

1.03262 

0.70598 

0.22943 

0.10783 

0.80894 

1.02319 

0.79060 

0.40426 

0.18999 

0.87770 

1.01496 

0.86477 

0.61701 

0.28998 

0.93282 

1.00826 

0.92517 

0.82273 

0.38666 

0.97229 

1.00342 

0.96897 

0.96378 

0.45295 

0.99470 

1.0006b 

0.99405 

1.00000 

0.46997 

1.00000 

1.00000 

l.OOOOO 

RE/lH-  1 

975318.1  RThI-  36147.1 

FRO-0. 920584. 

FC-1. 056867. 

<CF1-  2.07820. 

KCDI-  2.47275.  KCO-  2.33970.  TH-  0.04394 

Y-  6.90503. 

M-0. 269827.  DELTA- 

0.46997.  0-2-0-  0.06461. 

DELTA—  0.06343. 

DEL-A-  0.06372.  H-  1.443554.  N-  6.94305 

R-  6.96875. 

PO-  74.266.  PE- 

70.602.  TO-  474. 

929.  TAU-V- 

187.80  IN.. 

XSTAR-  .  37.561  IN..  HRH— 0.155787 

WAVE  END 

HbO 

.269827.  PO/POO-  0.99538.  PE/PEO-  0. 

99361.  TAO- 

563.41  IN.. 

DELAY  TIME-  0.094713  SECONDS.  V-  888.925 

WAVE  END 

6.94305.  H-  1.443554 

.  hRH— 0.155787. 

HSUM/VO-0. 

004499698.  HSUM/VE-0. 004499699.  WAVE  LGTH/TAO-l. 00000 

INVISCID 

ThROAT  m-0. 955942.  STREAMTUBE  THROAT  M-0 

.957060.  OEL-A/RaO. -0.000965.  OELTA/RAD. -0.059742.  TUBE-  24.66  FT. LG 

REYNOLDS  NO.(TAO)s  A.931&3D  OSt  OELTA/TaO-  B.3A1S3D-0A.  DEL*A/TAOs  1.13I020-OA.  NET  TIME-  O.OA20OO  SECONDS 


Y/OELTA 

Y(lNCH) 

U/UE 

T/TE 

MASS  FLUX 

0.0 

0.0 

0.0 

1.11178 

0.0 

0.00045 

0.00057 

0.35920 

1.07480 

0.33420 

0.00254 

0.00326 

0.45249 

1.06461 

0.42503 

0.01068 

0.01369 

0.54751 

1.05399 

0.51946 

0.03497 

0.04480 

0.64080 

1.04332 

0.61419 

0.09240 

0.11839 

0.72901 

1.03301 

0.70571 

0.20237 

0.25928 

0.80894 

1.02349 

0.79038 

0.37421 

0.47945 

0.87770 

1.01515 

0.86460 

0.59212 

0.75864 

0.93282 

1.00837 

0.92507 

0.80916 

1.03671 

0.97229 

1.00347 

0.96893 

0.96075 

1.23094 

0.99470 

1.00066 

0.99404 

1.00000 

1.28122 

1.00000 

1.00000 

1.00000 

RE/lNs  g97A25.t  RTHIs  9S738.f  FR0s0.919963»  FCsl.OSTAAS.  KCFIs  1.7990S.  KCDIs  2.11047.  KCD>  1.995821  TR-  0.10958 
APPR  5.  Y*  6.810711  Rsu. 2771581  DELTA-  1.2al22i  D-2-D-  0.16437i  DELTA—  0.15G25i  DEL*A-  0.15804i  H>  1.425919i  N-  7.53505 
R-  0.968751  PO-  74.2511  PE-  70.393T  TO-  474.903i  TAU^V-  539.48  IN.i  XSTAR-  109.808  IN.i  HRH— 0.155750 

WAVE  END  M-0. 2700321  PO/POO-  0.995181  PE/PEO-  0.99067i  TAO-  1006i09  IN.i  DELAY  TIME-  0.094713  SECONDS.  V-  888.760 
WAVE  END  N-  7.056361  H-  1.439320.  HRH—0. 155482.  HSUM/VO-0. 004482741 .  hSUM/VE-0. 004322653.  NAVE  LGTH/TAO-0. 69568 

INVISCIO  throat  M-0. 955942.  STREAMTUBE  THROAT  M-0. 958898.  DEL*A/RAO. -0.002565.  OELTA/RAO. -0.167056.  TUBE-  44.66  FT. LG 


REYNOLDS  NO. (TAO)-  9.028900  08.  OELTA/TAO-  1.2734TD-03.  DEL*A/TAO>  1.570860-04.  NET  TIME-  0.075000  SECONDS 
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NOMENCLATURE 

A*  Sonic  area 

Atube  Goss-section  area  of  a  Ludwieg  tube 
ao  Speed  of  sound  ahead  of  an  expansion  wave 

ai  Speed  of  sound  behind  an  expansion  wave 

C  Constant  in  Eq.  (29) 

Cf  Mean  skin-friction  coefficient 

Cf  Local  skin-friction  coefficient 

Fc  Ratio  of  incompressible  friction  coeffident  to  compressible  value 

FRs 

fin  Denotes  natural  logarithm  (base  e) 

log  Denotes  common  logarithm  (base  10) 

M  Mach  number 

n  Exponent  in  velodty  profile 

Pi  Static  pressure  behind  expansion  wave 

P4  Charge  pressure 

Re  Reynolds  number,  subscript  indicates  reference  length 

r  Radius  of  Ludwieg  tube^  also  recovery  factor  in  Eq.  (36) 

T  Temperature 

t  Time 

u  Velodty 

V  Velodty  of  concentrated  expansion  wave 

X  Distance,  Vt-x 
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X  Distance  along  Ludwieg  tube 

y  Distance  from  wall  of  Ludwieg  tube,  in. 

a  Variable  defined  by  Eq.  (38) 

P  Variable  defined  by  Eq.  (39) 

7  Ratio  of  specific  heats 

Ax  Length  traversed  by  head  of  expansion  wave  relative  to  end  of  Ludwieg  tube 

Aw  Length  from  head  to  tail  of  expansion  wave 

8  Boundary-layer  thickness 

5*  Displacement  thickness  of  boundary  layer  when  total  thickness  is  large  relative 

to  the  tube  radius 

di  Displacement  thickness  defined  by  Eq.  (9) 

5]c  Kinematic  displacement  thickness,  Eq.  (4S) 

r)  Nondimensional  variable  defined  by  Eq.  (5) 

6  Momentum  thickness  of  boundary  layer  when  total  thickness  is  large  relative  to 

tube  radius 

$1  Momentum  thickness  defined  by  Eq.  (10) 

dc  Equivalent  flat-plate  momentum  thickness 

dyi  Kinematic  momentum  thickness,  £q.  (42) 

K  Constant  in  logarithmic  velocity  profile 

M  Viscosity 

n  Wake  variable  in  logarithmic  velocity  profile 

p  Density 

p*  Density  integral  defined  by  Eq.  (13) 

Tw  Shear  stress  at  wall 
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SUBSCRIPTS 

aw  Adiabatic  wall  value 

i  Incompressible  value 

w  Wall  value 

1  Value  behind  expansion  wave  or  outside  of  boundary  layer,  except  5x  and  di 
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